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Abstract 

We investigate the cascading effects of extremely high energy (EHE) photons in 
the Earth's magnetosphere assuming that these photons arrive with the parameters of 
the highest energy AGASA events (energies, arrival directions). For the location of the 
AGASA Observatory, we determine the directions in the sky from which photons can 
cascade with a high (low) probability. In the case of the primary photons with the param- 
eters of the events with the energies > 10^° eV, we compute the average cascade spectra 
of secondary photons entering the Earth's atmosphere, and estimate their fluctuations 
around these average values by selecting the events with the largest and smallest number 
of secondary cascade photons. It is shown that most photons with the parameters of the 
highest energy AGASA events should initiate cascades in the Earth's magnetosphere with 
a high probability even though they tend to arrive from directions in the sky for which 
the perpendicular component of the magnetic field is weaker. On the other hand, if these 
events are caused by the photons with lower energies, then the fluctuations in their shower 
development in the magnetosphere and the atmosphere should be higher than in the case 
of photons with the energies estimated by the AGASA experiment. 

1 Introduction 

It is believed that the number of photons in the cosmic rays with the energies above 
10^^ eV may be significant. Extremely high energy cosmic ray hadrons should produce 
many photons in collisions with the Microwave Background Radiation when propagating 
through the universe (Wdowczyk & Wolfendale 1990, Halzen et al. 1990). The EHE 
photons may also be efficiently produced from the decay of massive particles (e.g. Higgs 
and gauge bosons), as predicted by some more exotic theories (see for review Bhattacharjee 
& Sigl 2000). Therefore, the investigation of the photon content in the EHE cosmic rays 
(CR) can give us important information about its origin. It is very difficult to obtain 
experimental constraints on the content of EHE CR due to the low statistics of detected 
particles at these energies and problems with the particle interaction models. In spite 
of these problems first efforts have been undertaken. For example, the analysis of the 



Haverah Park data allows the authors to put the upper limit on the photonic content in 
the cosmic rays above 4 x 10^^ eV equal to 55% (Ave et al. 1999, 2002). Similar limits have 
recently been obtained from the analysis of the AGASA data (Shinozaki et al. 2002). On 
the other hand, there is also experimental evidence, e.g. small scale clustering of EHE CR 
with the energies above 4 x 10^^ eV (Takeda et al. 1999, 2001, Tinyakov & Tkachev 2001), 
which suggest that the EHE cosmic rays contain a significant part of neutral particles (7- 
rays, neutrons ?). 

Even if photons may not completely dominate the highest energy cosmic rays, it is 
useful to investigate the hypothesis of their photonic origin since photons of sufficiently 
high energies may convert into pairs when propagating in a perpendicular magnetic 
field (Erber 1966). It has been shown by McBreen & Lambert (1981) that such a con- 
version of photons should also occur in the Earth's magnetic field. The secondary 
pairs created by photons should produce many secondary synchrotron photons, which in 
turn can again convert into pairs. Cascades initiated by the primary photons in the 
Earth's magnetosphere have already been discussed in the general context of the photon 
initiated cascade theory in the Earth's atmosphere by e.g. Aharonian, Kanevsky & Sa- 
hakian (1991), Karakula (1997), Kasahara (1997), Stanev & Vankov (1997), Anguelov & 
Vankov (2000) and also in the context of detecting high energy photons by planned Auger 
Observatories by Bednarek (1999) and Bertou, Billoir & Dagoret-Campagne (2000). 

Different experiments (Volcano Ranch, Haverah Park, Yakutsk, Fly's Eye, AGASA) 
have reported several showers with the energies > lO^^eV (Nagano & Watson 2000, Takeda 
et al. 1999, Hayashida et al. 2000). Best documented are the events observed by the 
AGASA array|^. This experiment registered 64 events with the energies > 4 x 10^^ eV 
(between them 10 have energies above 10^° eV) with zenith angles smaller than 45°. 
Extending the analysis of the AGASA data to the events coming for the zenith angles up 
to 60° increased the number of candidate events above 10^° eV to 17 (Sakaki et al. 2001). 
Unfortunately, the information which allow us to determine of the azimuth and zenith 
angles of these events, has only been published for 47 events with the energies > 4 x 10^^ 
eV, 8 of them have the energies > 10^° eV (Hayashida 2000). The parameters of two 
additional recent events with the energies above 10^° eV are reported on the AGASA web 
page. 

In this paper we analyse the features of cascades initiated by photons assuming 
that they arrive in the Earth's magnetosphere with the arrival directions of the AGASA 
events and energies > 5 x 10^^ eV. Our purpose is to find out if specific photons with the 
parameters of the AGASA events may interact with the Earth's magnetic field. 

2 Cascades initiated by photons in the Earth's mag- 
netic field 

The EHE photon with the energy can convert into an pair in the magnetic field B 
if the value of the parameter X7 = {E^/2mc^){B / Bcr) (where Bcr ~ 4.414 x lO^^G and 
mc^ is the electron rest mass) is high enough (Erber 1966). The secondary pairs can 
in turn produce synchrotron photons in the magnetic field. The efficiency of this process 
is determined by the parameter Xe± = {Ee/ 'rnc^){B / B^r) ■ The energies of synchrotron 

^available through AGASA web page: www-akeno.icrr.u-tokyo.ac.jp/AGASA/pub/ 



photons can be high enough to produce the next generation of e pairs, sustaining, in 
this way, the development of a cascade in the magnetic field. The simulations of the 
primary photons with energies of the highest energy cosmic ray events entering the Earth's 
magnetosphere show that up to a few generations of secondary pairs can be produced. 

We simulate the development of such a photon initiated cascade in the Earth's 
magnetic field by applying the spectra in the relativistic quantum domain for the pair 
production by photons and for the emission of synchrotron photons by secondary pairs 
given by Baring (1988). The spectra of pairs and synchroton photons computed using 
these formulae start to differ significantly in their high energy end in respect to the spectra 
without quantum corrections (e.g. Erber 1966) if the parameters X7 and Xe± sue much 
higher than one. In order to show the quantitative differences between the spectra with 
and without quantum corrections. Fig. |1] presents the probabilities of production of 
pairs by photons and the probabilities of emission of synchrotron photons by electrons 
(positrons) for selected values of x-y Xe±- These probabilities have been calculated 
from 

P = -^e,7 / Re,'rde, (1) 

where P is the random number, e is the ratio of a simulated energy of the electron 
(positron) normalized to the energy of the parent photon for the process of photon con- 
version into an pair, or the ratio of simulated energy of synchrotron photon normalized 
to the energy of the parent electron (positron) for the process of synchrotron emission of 
photon by an electron (positron). and Re are the spectra for production of pairs 
by the photon or for production of synchrotron photons by the electron (positron) (see 
Baring 1988). The mean free paths, Ae,7, for these processes are obtained from 

X-] = J\e,,de. (2) 

From Fig. |1], it is evident that if the values of x >> 1, the energies of produced secondary 
particles are significantly higher when the quantum corrections are included. Therefore, 
the more precise formulae given by Baring should be used if the parameters x-y ^i-nd Xe± 
are of the order of a few. This is the case of photons with the parameters of some highest 
energy AGASA events (> 10^° eV). With these exact formulae the cascade initiated by 
photons becomes more penetrating. 

To find out which specific photons can initiate cascades in the Earth's magnetic field, 
we have computed the mean free path for a conversion of photons with energies into 
pairs in a perpendicular magnetic field and for the production of synchrotron photons 
by electrons (positrons) with energies Ee in this same magnetic field. In Fig. |^ we show 
the results as a function of photon and electron energies in a perpendicular magnetic field 
with the strength of 0.3 G (typical of the Earth's surface on the equator). The mean free 
paths for other magnetic fields can be obtained from Fig. Q by simple linear scaling of 
the present computations (shift to the left and down for a stronger magnetic field). From 
Fig. 1^ it becomes clear that photons with the energies above ~ 3 x 10^^ eV may convert 
into pairs (the mean free path comparable to the Earth's radius). The shortest mean 
free paths (~ 45 km) have photons with the energies of a few 10^° eV. The mean free 
path for the production of synchrotron photon by an electron (positron) is equal to ~ 6 
km for Ee± < 10^^ eV. Our computations of the mean free paths for these two processes 
are consistent with the results obtained by Kasahara (1997). 



To analyse the cascade initiated by a specific photon with the parameters of the 
AGASA event we apply the Monte Carlo method. We start to follow the propagation 
of a photon from the distance of 5 Earth radii. The distance, x^, at which the photon 
converts into an pair is obtained from: 



where P2 is the random number. Note that depends on the propagation distance x. 
The energies of secondary electrons (positrons) are obtained from Eq. The distance, 
Xe±, of the emission of a synchroton photon by an electron (positron) is obtained from 
Eq. ^ in which A-y is replaced by the mean free path for emission of synchrotron photon Ae 
(see Eq. The energy of synchrotron photon produced by secondary electron (positron) 
is obtained from Eq. 0. 

Our simulations show that in the case of cascades initiated by photons with the 
energies of ~ 10^'' eV in the magnetic field with a typical value of 0.3 G, up to 4 pairs 
of secondary pairs can be produced. These 8 leptons produce on their path through 
the magnetosphere up to a few hundred synchrotron photons which arrive in the Earth's 
atmosphere with the energies above 10 GeV. Since leptons produced in the Earth's mag- 
netosphere seperate from each other (different charge) and their number is relatively low, 
the collective effects on emission of synchrotron radiation such as backreaction (Nelson & 
Wasserman 1991, Aloisio & Blasi 2002a) or coherent emission (Aloisio & Blasi 2002b) are 
neglected. 

We conclude that the photons with the parameters of the highest energy events 
observed in cosmic rays can convert into pairs in the Earth's magnetic field. However, 
due to the dipole structure of the Earth's magnetic field, the probability of conversion 
should strongly depend on the location of the cosmic ray observatory and on the arrival 
direction of the primary photon. 

3 Photons arriving to the location of the AGASA 
array. 

In order to have an idea which photons with the highest energy AGASA events can 
interact with the Earth's magnetosphere, we compute the profiles for the perpendicular 
components of the magnetic field for different photon arrival directions in respect to 
the location of the AGASA observatory (35°47'N, 138°30'E), using the World Magnetic 
Model0. The results of the calculations for different zenith Z and azimuth A angles 
(measured clockwise from the North) are shown in Figs. ^. It is evident that photons 
arriving to the AGASA array from the southern directions (the azimuths 90° < A < 
270°) meet the lowest perpendicular component of the magnetic field along their way 
of propagation. For these azimuth angles, and at some range of the zenith angles, the 
perpendicular component of the magnetic field can reach very small values since in these 
directions of the sky the magnetic field is almost tangent to the path of the photon (see 
characteristic cusps in Figs. |^). The magnetic field profiles determine the probability of 
conversion of a primary photon with specific energy into an pair. In Figs. ^ we show 

^available through: www.ngdc.noaa.gov 
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the probabilities for interaction of photons with different energies: = 2 x 10'^^ eV, 
10^° eV, and 5 x 10^^ eV, coming to the AGASA array from the sky. These probabihties 
strongly depend on the arrival directions and photon energies. Photons with the energies 
> 10^° eV, arriving from the northern sky for zenith angles greater than ~ 20°, should 
convert into pairs with the probability equal to one. However, for the southern sky, 
there is a region centered on Z ~ 45° and A ~ 175° (slightly shifted from the South 
to the East due to the south-east Asian magnetic anomaly), where the probability of 
interaction is significantly lower than one even for the highest energy AGASA events. 
The probability of photon interaction drops very fast with photon energy. For example, 
photons with the energies of 5 x 10^^ eV already have a chance to arrive in the Earth's 
atmosphere without interaction even from the northern sky. For such photons the region 
of low interaction probability expends largely containing in its center the region with the 
probability of interaction equal to zero. 

3.1 Cascading effects of photons with the parameters of the 
AGASA events. 

Let us assume that the highest energy AGASA events are caused by the EHE photons. 
For photons with the energies of some of these events the parameter Xy can be as high as 
a few. Therefore these photons should pre-cascade in the Earth's magnetosphere before 
developing the cascade in the Earth's atmosphere. Published parameters of the AGASA 
events with the energies above 5 x 10^^ eV (Hayashida et al. 2000) allows us to re- 
calculate their arrival directions, i.e. the zenith and azimuth angles (see Table |I]). Using 
them, we simulate cascades initiated by photons in the Earth's magnetosphere. The 
probabilities that photons will cascade with the parameters of these events (averaged 
over 10^ simulations) and the numbers of secondary photons with energies > 10^'' eV 
(averaged over 200 simulations) arriving at the distance of 20 km from the Earth's surface 
are reported in the fifth and the sixth column of Table |l|. First photon interaction usually 
occurs closer to the detector than the one Earth's radius. The most probable distance of 
the first interaction is usually closer than ~ 2500 km from the AGASA array (the results 
for some specific events are shown in Fig. |^). 

It is interesting to note that most of the AGASA events with the energies > 10^° 
eV (eight out of ten) tend to arrive from the southern part of the sky (their directions 
are marked by crosses in Figs. ^), i.e from the part of the sky where the probabilities of a 
photon conversion into an pair are the lowest. However, our simulations show that the 
photons with the parameters of these events cascade efficiently in the Earth's magnetic 
field. The probabilities of photon conversion into an pair are very high for most events 
with energies > lO^*' eV (see Table |1|). The photons with the parameters of three AGASA 
events, which have energies in the range (9 — 10) x 10^^ eV, seem to behave similarly to 
the photons with energies > 10^° eV. These two lower energy events, which arrive from 
the southern directions, also have a high probability of cascading. 

The situation with lower energy events, i.e. the events with the energies of (5 — 8) x 
10^^ eV, is different. Their arrival directions are quite isotropic in the azimuth since nine 
out of twenty events arrive from the southern sky and the rest from the northern sky. Only 
the events arriving from the northern sky cascade with a high probability (see Table |l|). 
For example, two AGASA events with the same energies 5.53 x 10^^ eV (91/05/31 and 
92/02/01), arriving from the northern and southern directions cascade with a completely 



different probability (equal to 0.935 and 0.005, respectively). Therefore, basing on the 
features of the cascades initiated by the photons with parameters of the AGASA events we 
distinguish two groups of particles. The first one contains events with the energies above 
9 X 10^^ eV (13 events) which tend to arrive from the southern directions. The second 
one contains events with the energies below 8 x 10^^ eV (20 events). Their distribution 
of arrival directions in the azimuth is more uniform. These two groups are separated by 
the lack of events with the energies in the range (8 — 9) x 10^^ eV. 

In Table ^we also report the number of secondary photons with the energies > 10^'' 
eV produced in cascades initiated by these events. These secondary photons earring most 
of the energy of a primary photon, fall into the Earth's atmosphere within a very small 
cone with the radius of a few tens of centimeters. They develop electromagnetic showers 
which in fact cannot be distinguished by cosmic ray detectors from a single shower. The 
secondary electrons and positrons produced in the cascade (usually not more than eight 
particles) arrive in the Earth's atmosphere with typical energies above ~ 10^^ eV. They 
carry only a few percent of the primary photon energy so their contribution to the cascade 
in the atmosphere is negligible. 

We also investigate the spectra of secondary cascade photons, produced in the mag- 
netic field by primary photons with the parameters of the AGASA events with the energies 
> 10^° eV, and fluctuations of these spectra in respect to the spectrum averaged over 200 
simulations. In Fig. ^ we show the histograms with the numbers of secondary photons 
within the specific energy range (A^ = AN^/ A{logE^), where A{logE^) = 0.1), which 
arrive in the atmosphere (i.e. at the distance of 20 km from the surface). The results of 
the simulations with the lowest and the largest numbers of secondary photons (selected 
from these 200 simulations) are marked by the dotted and dashed histograms. In the case 
of four AGASA events (with energies: 2.46 x 10^° eV, 2.13 x 10^° eV, 1.34 x 10^° eV, 
and 1.04 X 10^*^ eV) all 200 simulated photons converted into pairs. The fluctuations 
between different simulations in the number of secondary photons with energies > 10^^ 
eV may be as high as a factor of 2-3 around the average values. However, these fluctua- 
tions mainly concern the lower energy secondary photons and therefore, do not affect the 
development of cascades in the atmosphere seriously. The fluctuations in the location of 
the maximum in the spectrum are more important since this maximum is usually located 
close to or just above ~ 10^^ eV. Above this energy the Landau-Pomeranchuk-Migdal 
(LPM) (Landau & Pomeranchuk 1935, Migdal 1956) effect can have a significant influ- 
ence on their further cascading in the atmosphere. The fluctuations of the location of 
the maximum in the spectrum of secondary photons strongly depend on the considered 
AGASA event. They are only of the order of two in the case of events with the energies of 
2.13 X 10^° eV and 1.34 x 10^° eV. However, they can be as high factor of four for the 
event with the energy of 1.04 x 10^° eV. In fact, the primary photons with the parameters 
of some events with the energies > 10^° eV may not cascade at all in the magnetosphere. 
This possibility takes place rather rare, since all these events, except the event with the 
energy of 1.2 x 10^° eV, cascade with a probability greater than a half. 

A closer inspection of Fig |^ shows that most primary photons with parameters of 
the AGASA events with the energies > 10^° eV arrive in the atmosphere as a bunch 
of lower energy secondary photons with the mean energy just above 10^^ eV. Therefore, 
these primary photons should initiate cascades in the Earth's atmosphere which are not 
influenced by the LPM effect. The LPM effect is responsible for large fluctuations of 
showers which are induced by photons with the energies above a few 10^^ eV, as recently 



verified in the Monte Carlo simulations based on the Cosmos code (Kasahara 1997) and 
AIRES program (Bertou, Billoir & Dagoret-Campagne 2000). Also the depth of the 
maximum of showers simulated with the LPM effect is much deeper in the atmosphere 
than in the case of showers without the LPM effect. That is why it is generally expected 
that showers in the atmosphere initiated by the EHE photons with the energies above a 
few 10^^ eV should be easily distinguished from the ones induced by hadrons. However, 
as we noted above, primary photons with parameters of the AGASA events with energies 
above > 10^° eV usually pre-cascade in the magnetosphere with a high probability. Their 
products arrive in the atmosphere in the form of bunches of secondary photons initiating 
cascades in the atmosphere which are not significantly influenced by the LPM effect. 
The depth of the shower maximum which is induced by photons with the energies of 
~ 10^^ eV (equal to ~ 900 g cm^) does not differ significantly from the depth of shower 
maximum induced by the protons with energies of AGASA events (Bertou, Billoir & 
Dagoret-Campagne 2000). 

The photons with parameters of the lower energy AGASA events < 8 x 10^^ eV, 
which arrive from the southern directions, do not cascade in the Earth's magnetosphere 
or cascade with a low probability (Table |l|). The cascades initiated in the atmosphere 
by these primary photons, which do not suffer pre-cascading in the magnetosphere, are 
susceptible to large fluctuations due to the LPM effect. Hence, these showers can have 
the shower maximum even deeper in the atmosphere than showers initiated by primary 
photons with significantly higher energies which pre-cascaded in the magnetosphere. 

3.2 Cascading effects of photons with reduced energies of the 
AGASA events 

The conclusion mensioned in the previous section suggests that also the photons which 
have the arrival directions of AGASA events with energies > 10^° eV, but energies reduced 
by a factor of two or three in respect to the values estimated by AGASA group, should 
arrive in the atmosphere without pre-cascading in the magnetosphere. To check this we 
simulate the cascades initiated by primary photons with the parameters of the highest 
energy AGASA events but with energies reduced by factors mentioned above. It is found 
that in most cases the probability of conversion of such photons becomes significantly 
lower (only two events, 01/05/10 and 94/07/06, have the probability of interaction still 
greater than a half, see columns seven and eight in Table 

As an example, in Fig. we compare the average spectra of secondary photons 
produced by primary photons with reduced energies with the ones produced by primary 
photons with the energies reported by the AGASA experiment for two events: 93/12/03 
and 96/01/11. In both cases the average energies of secondary photons are higher in the 
case of primary photons with reduced energies than in the case of the average energies 
of secondary photons produced by primary photons with energies originally estimated by 
AGASA group (2.13 x 10^° eV and 1.44 x 10^° eV). Therefore we conclude that the highest 
energy AGASA events with energies reduced by a factor of two or three are more sensitive 
to fluctuations during cascading in the atmosphere because of a larger influence of the 
LPM effect on the development of showers initiated by the secondary photons with higher 
energies. We suggest that due to these fluctuations, not taken into account in derivation 
of the initial energies of the AGASA events, the energies of AGASA events arriving from 
the southern directions may be incorectly estimated. 



4 Conclusions 



We analyse the AGASA events with energies above 5 x 10^^ eV for which the arrival 
parameters have been published (Hayashida et al. 2000 and AGASA web page). It seems 
that two groups of events can be distinguished. The first group contains 13 events with 
energies > 9 x 10^^ eV. These events tend to arrive in the AGASA array mainly from the 
southern directions (10 out of 13). The second group contains 20 events with energies 
< 8 X 10^^ eV. Their arrival distribution in the azimuth angle is uniform (9 from the 
southern part of the sky and 11 from the northern sky). These two groups are separated 
by the lack of events with energies between (8 - 9) x 10^^ eV. 

We have shown that primary photons with the parameters of the highest energy 
AGASA events (> 9 x 10^^ eV) usually efficiently pre-cascade in the Earth's magneto- 
sphere in spite of the fact that most of them arrive from the southern directions where 
the perpendicular component of the magnetic field is the lowest. Photons with parame- 
ters of the lower energy AGASA events, (5 — 8) x 10^^ eV, can cascade efficiently when 
arriving from the northern directions. They usually do not cascade (or cascade with a 
low probability) when arriving from the southern directions. 

As we have mentioned the primary photons with the energies of 5 — 8 x 10^^ eV 
arriving from the southern direction should not cascade (or cascade with a low probabil- 
ity). They fall onto the atmosphere with average energies higher than the average energy 
of the secondary photons produced in cascades initiated in the Earth's magnetosphere by 
the primary photons with the energies > 9 x 10^^ eV. Therefore, the lower energy primary 
photons, (5 — 8) X 10^^ eV, should initiate showers in the atmosphere which fluctuate 
stronger due to the LPM effect than the showers initiated by the primary photons with 
the higher energies (> 9 x 10^^ eV). We suggest that these effects may influence the estima- 
tion of energies of photonic events arriving from the southern direction. This suggestion 
is consistent with the tendency that the higher energy AGASA events (> 9 x 10^^ eV) 
more likely arrive from the southern directions which is not the case in the lower energy 
events (< 8 x 10^^ eV). Their arrival distribution in the azimuth angle is consistent with 
the uniformity. 
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Table 1: The probabilities of photon interaction. 
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Figure 1: The probabilities of production of an pair with energy by the photon 
with energy (top figures) and for emission of synchrotron photons with energy E^gyn by 
the electron with energy E^ (bottom figures) for the selected values of the parameters x-y 
and Xe- The computations with the Baring's rates are shown by the full curves and with 
the Erber's rates by the dashed curves. 
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Figure 2: The dependence of the mean free paths for photon conversion with energy 
into an pair (full curve) and for production of the synchrotron photon by an electron 
(positron) with energy £'e± in the perpendicular magnetic field with the strength of 0.3 
G on energy of the primary particle. 
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Figure 3: The magnetic field profiles along the direction of EHE photons motion at the 
location of the AGASA Observatory (the distance measured in units of the radius of the 

Earth Rz). The left figure shows the magnetic field profiles as a function of the zenith 
angle: Z = 0° (full); A = 0°, and Z = 30° (dotted), and Z = 60° (dashed); and A = 180°, 
and Z = 30° (dot-dashed), and Z = 60° (dot-dot-dot-dashed), and the right figure for 
Z = 45° as a function of azimuth angle; A = 0° (full), 45° (dotted), 90° (dashed), 135° 
(dot-dashed), and 180° (dot-dot-dot-dashed). 
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Figure 4: The probability of conversion of photons with energies 2 x 10^° eV, 10^° eV, and 

5 X 10^*^ eV, arriving in the AG AS A array from different directions defined by the zenith 
Z and azimuth A angles. The probability, P, is less than one inside the region marked by 

fk^i A, 11 1^Q,-1^ P ^ n /.-looViQi-l^ P ^ n 1 ( Ar^\^\^i^A\ on,-1 
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Figure 5: The probability (AN/Alog(R/Rz)) of the first interaction of photons with 

parameters of some AGASA events (2.13 x 10^° eV - full histogram, 1.5 x 10^° eV - 
dotted, 1.2 X 10^° eV - dot-dashed, 1.04 x 10^° eV - dot-dot-dot-dashed, and 1.01 x 10^° 
eV - dashed) as a function of distance from the surface of the Earth, measured in units of 
the Earth's radius Rz- The histograms are averaged over 10^ simulated primary photons. 




Figure 6: Average energy distributions of secondary photons (full histograms) coming 

in the Earth's atmosphere. These secondary photons are produced in cascades initiated 

by the primary photons with the parameters of the AGASA events with the energies 

> 10^° eV (for other parameters sec Table 1). The dotted and dashed histograms show 

the simulations with the smallest and the largest numbers of secondary photons selected 
onn 




Figure 7: Average energy distributions of secondary photons (full histograms), produced 
by primary photons with parameters of the AGASA events 93/12/03, = 2.13 x lO^'' 
eV (on the left), and 96/01/11, = 1.44 x 10^° eV (on the right), averaged over 200 
simulated primary photons. The dashed and dotted histograms show the cases with the 
arrival directions of these events but energies reduced by a factor of 2 and 3, respectively. 



